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Abstract 



We develop and discuss an infrared-finite factorization and optimized renor- 
malization scheme for calculating exclusive processes which enables the inclu- 
sion of transverse degrees of freedom without entailing suppression of cal- 
culated observables, like form factors. This is achieved by employing an 
analytic, i.e., infrared stable, effective coupling a s (Q 2 ) which removes the 
Landau singularity at Q 2 = Aq CD by a power-behaved correction. The ensu- 
ing contributions to the cusp anomalous dimension, related to the Sudakov 
form factor, and to the quark anomalous dimension, which controls evolution, 
lead to enhancement of the hard part of exclusive amplitudes, calculated in 
perturbative QCD. The phenomenological implications of this framework are 
analyzed by applying it to the pion's electromagnetic form factor and the 
pion-photon transition. 

ll.10.ffi, 12.38.Bx, 12.38.Cy, 13.40.Hq, 13.40.Gp 
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1. Introduction 



In the last few years, several works ]^J3,|j,E| (among many others) addressed the pos- 
sibility of power corrections to the strong running coupling, beyond the operator product 
expansion. Such corrections, which are subleading in the ultraviolet (UV) region, correspond 
to nonanalytical contributions to the /3-function. 

The existence of these power corrections, if proven true, would greatly affect our under- 
standing of nonperturbative QCD effects. For instance, a power correction to a s gives rise 
to a linear term in the interquark static potential . On a more speculative level, one may 
argue Q that the source of such terms are small-size fluctuations in the nonperturbative 
QCD vacuum, perhaps related to magnetic monopoles in dual QCD or nonlocal condensates. 
Besides, and in practice, a power-behaved contribution at low scales can be used to remove 
the Landau singularity, present in perturbation theory, supplying in this way an infrared 
(IR) stable effective coupling . 

The aim of the present work is to develop in detail a factorization and renormalization 
scheme, which self-consistently incorporates such a nonperturbative power correction, and 
then use it to assess and explore exclusive processes. We do not, however, propose to involve 
ourselves in the discussion of whether or not such power corrections have a foundamental 
justification. We consider the ambiguity in removing the Landau pole as resembling the 
ambiguity in adopting a particular (non-IR-finite) renormalization scheme in perturbative 
QCD. The justification for such a procedure will be supplied a posteriori by the self-consistent 
incorporation of higher-order perturbative corrections and by removing any IR- sensitivity 
of calculated hadronic observables. 

A key ingredient of our approach is that the modified effective coupling will be taken 
into account not only in the factorized short- distance part, i.e., through the fixed-order 
perturbation expansion, but also in the resummed perturbative expression for soft-gluon 
emission and in the renormalization-group controlled evolution of the factorized parts. 

To this end, we adopt as a concrete power-corrected effective coupling, an analytic model 
for a s , proposed by Shirkov and Solovtsov ||, which yields an IR-finite and universal effec- 
tive (running) coupling. This model combines Lehmann analyticity with the renormalization 
group to remove the Landau singularity at Q 2 = Aq CD , without employing adjustable pa- 
rameters, just by modifying the logarithmic behavior of a s by a simple (nonperturbative) 
power term (minimality of the model). 

At the present stage of evidence, it would be, however, premature to exclude other 
parametrizations, since neither the sign of the power correction nor the size of its exponent 
are established, and one could introduce further modifications ||,[7j . 

Continuing our previous exploratory study || , we further extend and test our theoretical 
framework by also including into the calculation of the pion form factor the next-to-leading 
order (NLO) perturbative contribution to the hard scattering amplitude (see, e.g., and 
earlier references cited therein). 

The ultimate goal of the present analysis is not to obtain results in perfect agreement 
with the data, but to expose and discuss the conceptual advantages of our scheme relative 
to previous conventional approaches HlQjll,|l2| . 

The major advantage of such a theoretical framework, the object of this paper, is that 
it enables the inclusion of transverse degrees of freedom, primordial (i.e., intrinsic) 0] and 



3 



those originating from (soft) gluonic radiative corrections [10|, without entailing suppression 
of perturbatively calculated observables, viz., the pion form factor. This enhancement is due 
to power-term generated contributions to the anomalous dimensions of the cusped Wilson 
line, related to the Sudakov form factor, and such to the quark wave function which governs 
evolution. 

Although most of our considerations refer to the pion as a case study for the proposed 
framework, the reasoning can be extended to describe three-quark systems as well. This will 
be reported elsewhere. 

The outline of the paper is as follows. In the next section we briefly discuss the essential 
features of the IR-fmite effective coupling. In sect. [3] we develop and present our theoretical 
scheme. Sect. |] extends the method to the next-to-leading order contribution to the hard- 
scattering amplitude. In sect. [| we discuss the numerical analysis of the electromagnetic 
pion form factor, revolving around the appropriate kinematic cuts in the evaluation of the 
IR-modified Sudakov form factor which comprises additional nonleading contributions. We 
also provide arguments for the appropriate choice of the renormalization scale. In this section 
we also show the theoretical prediction for the pion-gamma transition form factor, derived 
with our theoretical framework, as an independent justification of the approach. Finally, in 
sect. ^ we summarize our results and draw our conclusions. 



2. Model for QCD effective coupling 

The key element of the analytic approach of Shirkov and Solovtsov is that it combines a 
dispersion-relation approach, based on local duality, with the renormalization group (RG) 
to bridge the regions of small and large momenta, providing universality at low scales. The 
approach is an extension to QCD of a method originally formulated for QED ||13|| . 

At the one-loop level, the ghost singularity is removed by a simple power correction and 
the effective coupling reads 



ln(Q 2 /A 2 ) A 2 - Q 2 \ ' w 

where A = Aqcd is the QCD scale parameter. 

This model has the following interesting properties. It provides a nonperturbative reg- 
ularization at low scales and leads to a universal value of the coupling constant at zero 
momentum a^\Q 2 = 0) = 47r//3 ~ 1.396 (for three flavors), defined only by group con- 
stants. No adjustable parameters are involved and no implicit "freezing", i.e., no saturation 
hypothesis of the coupling constant is invoked. 

Note that this limiting value (i) does not depend on the scale parameter A - this being 
a consequence of RG invariance - and (ii) extends to the two-loop order, i.e., af'(Q 2 = 
0) = a^(Q 2 = 0) = a s (Q 2 = 0). (In the following the bar is dropped.) Hence, in 
contrast to standard perturbation theory, the IR limit of the coupling constant is stable, 
i.e., does not depend on higher-order corrections and is therefore universal. As a result, 
the running coupling constant also shows IR stability. This is tightly connected to the 
nonperturbative contribution oc exp(— 47r/a/3o) which ensures analytic behavior in the IR 
domain by eliminating the ghost pole at Q 2 = A 2 , and extends to higher loop orders. Besides, 



47T 



^(Q 2 ) = f- 
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the stability in the UV domain is not changed relative to the conventional approach, and 
UV perturbation theory is preserved. 

At very low-momentum values, say, below 1 GeV, Aqcd in this model deviates from 
that used in minimal subtraction schemes. However, since we are primarily interested in a 
region of momenta which is much larger than this scale, the role of this renormalization- 
scheme dependence is only marginal. In our investigation we use Aq 11 ^ -3 '' = 242 MeV which 

corresponds to Aqcr? f=3) = 200 MeV. 

The extension of the model to two-loop level is possible, though the corresponding ex- 
pression is too complicated to be given explicitly 0. An approximated formula with an 
inaccuracy less than 0.5% in the region 2.5 A < Q < 3.5 A, and practically coinciding with 
the exact result for larger values of momenta, is provided by [0] 



a. 



(2) 
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(2) 



where = 11 - f n f = 9, fa = 102 - f n f = 64, and D x = 0.035 for n f = 3. 

With experimental data at relatively low momentum-transfer values for most exclu- 
sive processes, reliable theoretical predictions based on perturbation theory are difficult to 
obtain. Both the unphysical Landau pole of a s and IR instability of the factorized short- 
distance part are affecting such calculations. It is precisely for these two reasons that the 
Shirkov-Solovtsov analytic approach to the QCD effective coupling can be profitably used 
for computing amplitudes describing exclusive processes, like form factors, [|14],|15|,|nj . The 
improvements are then: (i) First and foremost, the nonperturbatively generated power cor- 

via the cusp anomalous dimension 



rection modifies the Sudakov form factor |p^7| , p!8| , p!9| , |20| , pT| 
22 1, and changes also the evolution behavior of the soft and hard parts through the modified 
anomalous dimension of the quark wave function. This additional contribution to the cusp 
anomalous dimension is the source of the observed IR enhancement and helps taking into ac- 
count nonperturbative corrections in the perturbative domain, thus improving the quality of 
the predictions, (ii) Factorization is ensured without invoking the additional assumption of 
"freezing" the coupling strength in the IR regime by introducing, for example, an (external) 
effective gluon mass to saturate color forces at large distances, (iii) The Sudakov form factor 
does not have to serve as an IR protector against a s singularities. Hence the extra constraint 
of using the maximum between the longitudinal and the transverse scale, as argument of 
a s , proposed in |10| and used in subsequent works, becomes superfluous, (iv) The factoriza- 
tion and renormalization scheme we propose on that basis enables the optimization of the 



(arbitrary) constants which define the factorization and renormalization scales [17 , [18 , 23 , 23 
This becomes important when including higher-order perturbative corrections. 



3. Infrared-finite factorization and renormalization scheme 

Application of perturbative QCD is based on factorization, i.e., how a short-distance part 
can be isolated from the large- distance physics. But in order that observables calculated 
with perturbation theory are reliable, one must deal with basic problems, like the resumma- 
tion of "soft" logarithms, IR sensitivity, and the factorization and renormalization scheme 
dependence. 
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It is one of the purposes of the present work to give a general and thorough investigation 
of such questions. 

The object of our study is the electromagnetic pion's form factor in the spacelike region, 
which can be expressed as the overlap of the corresponding wave functions between the 
initial and final pion states: El 



(3) 



where we have assumed dominance of the valence qq state, with e q denoting the charge of 
the struck quark, and where 



1 



kj 
k 



(1 - x)q ± , 



struck quark 
spectators . 



(4) 



The wave function ipn{x, kjj is the amplitude for finding a parton in the valence state with 
longitudinal momentum fraction x and transverse momentum k^. 

In order to apply a hard-scattering analysis, we dissect the pion wave function into a soft 
and a hard part with respect to a factorization scale //p, separating the perturbative from 
the nonperturbative regime, and write (in the light-cone gauge A + = 0) 



fa (x, kj = C° ft k ± ) 9 (/4 - ki) + ^ ard (x, kj 9 (k] 



i4 



(5) 



Then the large k± tail can be extracted from the soft wave function via a single-gluon 
exchange kernel, encoded in the hard scattering amplitude Th, so that [14] 

ip haxd (x, k ± ) = J dy J d 2 \ ± T R (x,y, ll;a a (lj.)) ip so ft(y, l±) ■ (6) 



As a result, the pion form factor in Eq. (|| 

^(Q 2 )=«®Cft + « 



is expressed in the factorized form 



Th ® Cft 



+ 



soft 



+ 



Cft ® Tb 



soft 



+ . . . 



(7) 



where the symbol <g> denotes convolution defined by Eq. (^|) . The first term in this expansion 
is the soft contribution to the form factor, with support in the low-momentum domain, that 
is not computable with perturbative methods. The second term represents the leading order 
(LO) hard contribution due to one-gluon exchange, whereas the last one gives the NLO 
correction. We will not attempt to calculate the first term here, but adopt for simplicity the 
result obtained by Kroll and coworkers in fl26|. For other, more sophisticated, attempts to 



model the soft contribution to F-^^Q 2 ), we refer to [p7| , p^1 . 

We now employ a modified factorization prescription [ p~Q] , p~T| , which explicitly retains 
transverse degrees of freedom, and define (see for illustration fig. [1]) 



hard 
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< ki < CU 2 Q 
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(8) 



with b, the variable conjugate to k±, being the transverse distance between the quark and the 
antiquark in the pion valence Fock state. The Sudakov-type form factor exp(— S) comprises 
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Figure 1. Illustration of the factorized pion form factor, exhibiting the different regimes of dynamics. 
The wiggly line denotes the off-shell photon. Gluon exchanges are not explicitly displayed. The region of 
hard-gluon rescattering (LO and NLO) is contained in the short-distance part, termed Tjj. The blobs e(— S) 
represent in axial gauge Sudakov-type contributions, whereas nonperturbative effects are absorbed into the 
pion wave functions ip™ and . 



leading and next-to-leading logarithmic corrections, arising from soft and collinear gluons, 
and resums all large logarithms in the region where <C Q 2 [[23|,^4|,[||] . The presence of 



these logarithms results from the incomplete cancellation between soft-gluon bremsstrahlung 
and radiative corrections. It goes without saying that the function S includes anomalous- 
dimension contributions to match the change in the running coupling in a commensurate 
way with the changes of the renormalization scale (see below for more details). 
Going over to the transverse configuration space, the pion form factor reads [|l(| 



dxdy 



d 2 b. 



o ./-oo (An) 2 



Vr (y, 6, P'; C lt C 2 , C 4 ) T H (x, y, 6, Q; C 3 , C 4 ) 



x V^{x,b,P;C u C 2 ,C A ) , 
where the modified pion wave function is defined in terms of matrix elements, viz., 

-|k x |<At 



(9) 



V n (x, b, P, fi) 



dz~ 
~2n~ 



d 2 k ]e - ik ^V w (x, k ± ,P) 
T (g(0) 7 + 7 5 g (o 



cP+z 







, z ,bj 



vr(P) 



A+=0 



(10) 



with P + = Qj\f2 = P~' i Q 2 = —{P' — P) 2 , whereas the dependence on the renormalization 
scale /i on the rhs of Eq. (|T0| ) enters through the normalization scale of the current operator 
evaluated on the light cone. (Note that we set all light quark masses equal to zero and work 
in the chiral limit, i.e., M n = 0.) 

A few comments on the scales involved: 

• The scale C3/6 serves to separate perturbative from nonperturbative transverse dis- 
tances (lower factorization scale). We assume that there exist some characteristic scale 
&nonp — (k5_) 1 / 2 /x(l — x) ~ 0.5 GeV, related to the typical virtuality (offshellness) of 
vacuum quarks. This scale should also provide the natural starting point for the evo- 
lution of the pion wave function. In the following, we match the nonperturbative scale 
C 3 /b with the scale Ci/b, where the resummation of soft gluons starts, i.e., we set 
C\ = C 3 . The lower boundary of the scale C\jb is set by Aq CD , though the results are 
not very sensitive to using a somewhat larger momentum scale, as we shall see later. 
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• The resummation range in the Sudakov form factor is limited from above by the scale 
C^Q (upper factorization scale). (Note that the constant here differs in notation 
by a factor of y/2 relative to that used by Collins, Soper, and Sterman [p^j , i.e., 
Cv? ss = v^2C2-) This scale may be thought of as being an UV-cutoff for the (soft) 
Sudakov form factor, and enables this way a RG-controlled scale dependence governed 
by appropriate anomalous dimensions within this subsector of the full theory. 

• Analogously to these factorization scales, characterized by the constants C\,C2, and 
C3, we have introduced an additional arbitrary constant C4 to define the renormaliza- 
tion scale Cif(x,y)Q = /ir, which appears in the argument of the effective coupling 
af 1 . The effective coupling plays a dual role: it describes the strength of the interac- 
tion at short distances, and controls via the anomalous quark dimension the lower and 
upper boundaries, respectively, for the evolution of Th, and V v to the renormalization 
scale. 

The appropriate choice of the unphysical and arbitrary constants Cj will be discussed in our 
numerical analysis in sect. |5|. 

The ambiguities parametrized by the scheme constants Cj emerge from the truncation 
of the perturbative series and would be absent if one would be able to derive all-order 
expressions in the coupling constant. In fact, the calculated (pion) form factor depends 
implicitly on both scales: the adopted renormalization scale via a s , and the particular 
factorization scheme through the anomalous dimensions. Since the latter also depend on a s , 
the factorization-scheme and the renormalization-scheme dependences are correlated. On 
the other hand, the physical form factor is independent of such artificial scales, and satisfies 

fi — — - = 0, for /i being any internal scale. Obviously, both scheme dependences should 
be treated simultaneously and be minimized in order to improve the self-consistency of 
the perturbative treatment. In order to render the perturbative prediction reliable, the 
parameters C{ should be adjusted in such a way, as to minimize the influence of higher-order 
corrections, thus resolving the scheme ambiguity. However, in the present investigation we 
are not going to explicitly match the fixed-order NLO contributions with the corresponding 
terms in the resummed expression for the "soft" logarithms. This check will be included in 
future work. 

In Eq. (HI), Th is the amplitude for a quark and an antiquark to scatter collinearly via a 
series of hard-gluon exchanges, and, in LO perturbative expansion, it is given by 

T R (x,y,b,Q ]f i R ) = 8C F aT(^ R )K (^bQ) . (11) 

This result is related to the more familiar momentum-space expression 

rp f , , n \ 167tC f « s (^r) no x 
T R {x,y,k ± ,l ± ,Q,fi R ) = — — — ——2 (12) 

xyQ 2 + (k ± + ljj 

via the Fourier transformation 

/oo 
d 2 b ± T H (x, y, b, Q, /j, R ) exp [zb± • (kj_ + l±)} , (13) 
-00 

where use of the symmetry of ip^ under x <->• 1 — x = x has been made, and where C-p = 
(iV c 2 - l)/2iV c = 4/3 for SU(3). 
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The amplitude 



V n (x, b,P ~Q, C 1 , C 2 , n) = exp 



-s (x, b, Q, d, C 2 ) - s (x, b, Q, d, C 2 ) 

-2 r ^j q «m 

JC\/b /i 



V v {x,b,C x /b) (14) 
describes the distribution of longitudinal momentum fractions of the qq pair, taking into 



account the intrinsic transverse size of the pion state JTTj] , and comprising corrections due 
to soft real and virtual gluons ||10|| , including evolution to the renormalization point. 

The pion distribution amplitude evaluated at the factorization scale is approximately 
given by 



V^x^Cr/b)-^ (x,d/6)S(x,6) . 



(15) 



In the present work, we follow Jakob and Kroll |T]J and parametrize the distribution in the 
intrinsic transverse momentum (or equivalently the interquark transverse distance b) in 
the form of a nonfactorizing in the variables x and k^ (or x and b) Gaussian function which 
is normalized to unity, 



S(x, b) = An exp 



Ag{xW 



(16) 



where g(x) = 1/xx, with an appropriate width (3 to be specified below. 

Neglecting transverse momenta in Eq. ( |T2D (collinear approximation), the only depen- 
dence on kj_ resides in the wave function. Limiting the maximum value of k±, these degrees 
of freedom can be integrated out independently for the initial and final pion states to give 
way to the corresponding pion distribution amplitudes, which depend only implicitly on the 
cutoff momentum: 



<f> n (x, n' 



m 2 d 2 k i 



167T 3 



(17) 



where /„- = 130.7 MeV and iV c = 3. Integrating on both sides of this equation over x 



normalizes (j)^ to unity, i.e., Jq 1 dx 4> T (x, /x 2 
leptonic decay 7r — > /U + z^ for any factorization scale. 



1 because the rhs is fixed to 



JjL 



by the 



Hence, the full (model) wave function for the pion takes finally the form |TT 



^ (x, kj 



<P(x)f3' 2 g(x)exp -g(x)f3 2 k 2 



(18) 



Let us now return to Eq. ( filf) . The Sudakov form factor F$ (£, b, Q, C\, C 2 ), i.e., the 
exponential factor in front of the wave function, will be expressed as the expectation value 
of an open Wilson line along a contour of finite extent, C, which follows the bent quark line 
in the hard-scattering process from the segment with direction P to that with direction P' 
after being abruptly derailed by the hard interaction which creates a "cusp" in C, and is to 
be evaluated within the range of momenta termed "soft" , confined within the range limited 
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by C\/b and (where £ = x, x, y, y). (This means that the region of hard interaction of 



the Wilson line with the off-shell photon is factorized out.) Thus we have p7|Jl^p0pi| , p0 



F s (W(C)) = (P exp Ug dz ■ t a A a (z) 



soft 



(19) 



where P denotes path ordering along the integration contour, and where < ... >a denotes 
functional averaging in the gauge field sector with whatever this entails (ghosts, gauge 
choice prescription, Dirac determinant, etc.). Having isolated a subsector of the full theory, 
where only gluons with virtualities between Ci/b and C 2 iQ are active degrees of freedom, 
quark propagation and gluon emission can be described by eikonal techniques, using either 



Feynman diagrams p^fl8|1 or by employing a worldline casting of QCD which reverts the 
fermion functional integral into a path integral 



Then the Sudakov functions, entering Eq. (13), can be expressed in terms of the 
momentum-dependent cusp anomalous dimension of the bent contour to read 



1 f c ^Q da 

s (£, b, Q, d, c 2 ) = - r cusp ( 7 , <») 

£ J Oi lb IX 



c\ I - cusp 

2 Jd/b Li 

with the anomalous dimension of the cusp given by 



(20) 



r cusp ( 7 ,<») 



(21) 



7 = In ( j being the cusp angle, i.e., the emission angle of a soft gluon and the bent 
quark line after the external (large) momentum Q has been injected at the cusp point by 
the off-mass-shell photon. The functions A and B are known at two-loop order: 



lK {al^))+{5{g)^K{C x ,aT{^)) 



c F <W) + W)c F 

ix 2 



n' 



7T / 



(22) 



and 



B «(M)) = "2 I* (Ci, O)) + G (£, C a> a^))} 



3 7T "'\Ci 4 



(23) 



respectively. The K-factor in the MS scheme to two-loop order is given by |I^JT9|j23]j2^j31 
K (&) = {% ~ 4- J C A - ^n f T F + & In (de^/2) (24) 



with C*a = A^ c = 3, rif = 3, T F = 1/2, j E being the Euler-Mascheroni constant. 

The quantities /C, Q are calculable using the nonabelian extension to QCD [23] of the 
Grammer-Yennie method |B2| for QED. Alternatively, one can calculate the cusp anomalous 
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dimension employing Wilson lines 0,^,^T|,Q. In this latter approach (see, e.g., pl|), the 
IR behavior of the cusped Wilson line is expressed in terms of an effective fermion vertex 
function whose variance with the momentum scale is governed by the anomalous dimension 
of the cusp within the isolated effective subsector. Since this scale dependence is strictly 
restricted within the low-energy sector of the full theory, IR scales are locally coupled and the 
soft form factor depends only on the cusp angle which varies with the interquark transverse 
distance b between C\/b and C 2 iQ- 

The corresponding anomalous dimensions, are linked to each other (for a nice discus- 
sion, see, |T7[) through the relation 2r cusp « n (yu)) = 7^ « n (yu)) with r cusp (a^ n (/i)) = 
C F a^ n (fi 2 )/7i, which shows that = A (af^/x)). (Note that j g = —7^.) 

The soft amplitude V n (x, b, Ci/b, /x) and the hard-scattering amplitude Th (x, y, b, Q, /x) 
satisfy independent RG equations to account for the dynamical factorization (recall that 
both b and £ are integration variables) with solutions controlled by the modified "evolution 
time" (see, e.g., and earlier references cited therein): 



T 



d 



dk 2 a^(k 2 



c 2 /b 2 k 2 



47T 



1 ^ In (/x 2 /A 2 ) 
A) n In (C\l (6A) 2 



1 

Jo 



In 



(Ci/by 



-In 



A 2 



91 

b 2 



A 2 



(25) 



from the factorization scale C\jb to the observation scale /x, with A denoting Aqcd as before. 
The evolution time is directly related to the quark anomalous dimension, viz., 7 q (a^ n (fi)) = 
—a^ n (fi 2 )/ir. One appreciates that the second term in (p5| ) stems from the power-generated 
correction of the effective coupling and is absent in the conventional approach. At moderate 
values of /x 2 it is "slowing down" the rate of evolution. 

The leading contribution to the IR- modified Sudakov functions s (£, b, Q, d, d) (where 
£ = x, x, y, y) is obtained by expanding the functions A and B in a power series in af 1 , and 

( a an \ n / \ m 
In \j££,bQ\ , which can be transformed back 

into large logarithms In (Q 2 /k_|_ 2 ) in transverse momentum space. Employing equations ([!]) 
and (0), the leading contribution results from the expression 



s(^b,Q,C u C 2 



a 



an(2) 



1 /°*«M 21n (c*ff 

2 Jd/b /X [ \ /X J [ Ti 

+ <^M B (i) {Cl ,c 2 ) + o 

IT \ 71 



00 



a 



an(l) 



7T 



^ 2 > (d) 



(26) 



where Eq. (fj) is to be used in front of A^ l \ whereas the other two terms are to be evaluated 
with Eq. ([!]). The specific values of the coefficients B^' are 



A® (d) 
(d,d) 



d 

^ F x(d; 



2 p 27b-1" 



In 



a 2 



(27) 
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As now the power-correction term in a s gives rise to polylogarithms, a formal analytic 
expression for the Sudakov form factor is too complicated for being presented. In the 
calculations to follow, Eq. (|26| ) is evaluated numerically with appropriate kinematic bounds. 
Note that, neglecting the power-generated logarithms, we obtain an equation for the ordinary 
Sudakov function, which we write as an expansion in inverse powers of /3q to read 



s{t,b,Q,C 1 ,C 2 ) = - 



1 

0o 



(2A^Q + B^)\ni-2A^ (Q -b 



P 2 o v 
where the abbreviations QlOj| Q i 



b b , 
Q Q-b 



In 2b 



+ 



ln 2 (2g) - In 2 (26)] I , (28) 



In and 6 



In ^ have been used. 



This quantity differs from the original result given by Li and Sterman in [10|, and, though 



it almost coincides numerically with the formula derived by J. Bolz f33 |, it differs from that 
algebraically. 

All told, the final expression for the electromagnetic pion form factor at leading pertur- 
bative order in Th and next-to-leading logarithmic order in the Sudakov form factor has the 
form 



2 r°° 
F«{Q 2 ) = -nC ¥ fl dx dy 6d6a s an(1) W0as(x)0as(2/)exp 
o Jo Jo Jo 



6 2 (xx + yy) 



xK {^/xyQb) exp \-S (x, y, 6, Q, Ci, C 2 , C 4 )] 



(29) 



where 



S (x, y, 6, Q, Ci, C 2 , C 4 ) = s (x, 6, Q, C u C 2 ) + s (x, 6, Q, C u C 2 ) + (x <-> y) - 8r ^ 

(30) 

with t (Sl, /i R ) given by Eq. (||) and /i R = C 4 f(x, y)Q. 

Before we go beyond the leading order in the perturbative expansion of the hard- 
scattering amplitude, Th, let us pause for a moment to comment on the pion wave function. 
We have proactively indicated in Eq. ( |2"9"D that its asymptotic form as (a;) = 6xx will be 
used. 

A few words about this choice are now in order. 

Hadron wave functions are clearly the essential variables needed to model and describe 
the properties of an intact hadron. In the past, most attempts to improve the theoretical 
predictions for the hard contribution to the pion form factor have consisted of using end-point 
concentrated wave functions. In this analysis, we refrain from using such wave functions 
of the Chernyak-Zhitnitsky (CZ) type |[34| , referring for a compilation of objections and 
references to [§] (see also |35|), and present instead evidence for an alternative source of 
enhancement due to the nonperturbative power correction in the effective coupling. This 
effect is found || to be quite significant, even for the asymptotic pion wave function which 
has its maximum at x — 1/2. Indeed, the IR-enhanced hard contribution can account 
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already at leading perturbative order for a sizable part of the measured magnitude of the 
electromagnetic pion form factor, though agreement with the data calls for the inclusion 
of the soft contribution (cf. Eq. (|7|)) P5,57,55|,26| - even if the NLO correction is taken 
into account (see sect. |5[). Nevertheless, the true pion distribution amplitude may well 
be a "hybrid" of the type $£ ue = 90%$^ s + 9%$^ z + \%cf /2) , where the mixing ensures 
a broader shape with the fourth-order, "Mexican hat" -like, Gegenbauer polynomial G\ , 
being added in order to cancel the dip of 0cz at x = 1/2. 

First tasks from instanton-based approaches show that the extracted pion distribution 
amplitudes are very close to the asymptotic form p9| , ^0| . Similar results are also obtained 
using nonlocal condensates fl4T|. The discussion of non-asymptotic pion distribution ampli- 



tudes will be conducted in a separate publication. 



4. Pion form factor to order (a^iQ 2 ))' 



Next, we generalize our calculation of the hard contribution to the pion form factor by 
taking into account the perturbative correction to Th of order a 2 , using the results obtained 
in [42,43,44,9|, in combination with our analytical scheme. 

To be precise, we only include the NLO corrections to Th leaving NLO corrections to the 
evolution of the pion distribution amplitude aside. The reason is that for the asymptotic 
distribution amplitude, at issue here, these corrections are tiny p3|fl. For subasymptotic 
distribution amplitudes, however, NLO evolutional corrections have to be taken into 
account. The calculation below does not claim to be strictly correct, the reason being that 
the transverse degrees of freedom in the NLO terms to Th have been neglected, albeit the 
intrinsic ones in the wave functions are taken into account. Hence our prediction should 
be regarded rather as an upper limit for the size of the hard contribution to the pion form 
factor than as an exact result. Taking into account the /cj_-dependence of Th at NLO, this 
result will be reduced, though we expect that due to IR enhancement this reduction should 
be rather small. Thus, we have 



16ttCV 



XK (y/xyQb) 0asO)0a S (2/) 

b 2 (xx + yy) 



x exp 



exp (-S (x, y, b, Q, C x , C 2 , C 4 )) 



x 



1 + — (iuv (x, <2 2 //4) + /ir (x,y,Q 2 /i4) + fc{x,y) 



(31) 



where the Sudakov form factor, including evolution, is given by Eq. 
the functions /j are taken from |9| 0. They are given by 



/ip = Ci/b, and 



We wish to thank Dr. B. Melic for sending us the corrected version of this paper prior to 
publication. 
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Ajv (x,y,Q 2 /n 
fm (x,y,Q 2 /fi 



(5 
4 \3 



- In (xy) + In ^ 
^ (3 + In (xy)) ( ~ In (xy) - In ^ 



/c (x, y) = — [-34 + 12 In (xy) + In x In y 
+ In x In y — In x In y — In x In y 
+(1 - x - y)H(x, y) + R(x, y)\ 



(32) 



and are related to UV and IR poles, as indicated by corresponding subscripts, that have 
been removed by dimensional regularization along with the associated constants ln(47r) — 7e, 
whereas the last term is scale-independent. In evaluating expression fc in fl32|), we found it 
particularly convenient to use the representation of the function H(x,y) given by Braaten 
and Tse 



H(x,y) 



1 — x — y 



[Li 2 (x) + Li 2 (y) — Li 2 (x) — Li 2 (y ) + In x In y — In x In y] 



(33) 



where Li 2 denotes the Spence function. For the function R(x, y) we have used the expression 
derived by Field et al. [j4"2"| , except at point x ~ y, where we employed the Taylor expansion 
displayed below: 



R(x,y) 



3(-l + y) y 2 



x-y 



3(-l+y) 2 y 3 
(x - y f 

30 (-1 + ^)3^4 



(-l + 33y-45y 2 + 13y 3 ) lny 
+ y (-l+y + (9-13y)ylny) 
~(-l + y) 2 (-l + lQy) lny 
+ y f-l + 13y-12y 2 + 2y 2 lny 



(-l+y) 6 (9-U8y + 9y 2 ) \ny 



(34) 



- y (9 - 148 y + 328 y 2 - 189 y 3 + y 3 (5 + 9 y) In ?/ 

Note that this expression does not reproduce its counterpart in |f42 . 

Having developed in detail the theoretical apparatus, let us now turn to the concrete 
calculation of the pion form factor up to NLO. 



5. Validity of the analysis 

This numerical analysis updates and generalizes our previous investigation in ||. In 
order to set up a reliable algorithm for the numerical evaluation of the expressions presented 
above, we have to ensure that this is done in kinematic regions where use of fixed-order 
or resummed perturbation theory is legal. Further, restrictions have to be imposed to 
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avoid double counting of gluon corrections by carefully defining the validity domain of each 
contribution to the pion form factor. These kinematic constraints are compiled below. 



Kinematic cuts 

1. Ci/b > Aq CD ; otherwise the whole Sudakov exponent exp(— S) (cf. Eq. fl30D ) is 
continued to zero because this large-6 region is properly taken into account in the 
wave functions. This condition excludes from the resummed perturbation theory soft 
gluons with wavelengths larger than C±/A, which should be treated nonperturbatively. 

2. C-£Q > Ci/b] otherwise each Sudakov exponent exp [— s (£, b, Q, C\, C2)] in Eq. (|3~0D is 
"frozen" to unity because this small-6 region is dominated by low orders of perturbation 
theory rather than by the resummed perturbation series, and consequently contribu- 
tions in this region should be ascribed to higher-order corrections to Th- Yet evolution 
is taken into account to match the scales in our "gliding" factorization scheme. 



3. C i f{x,y)Q > C\jb\ otherwise evolution time t (Ci/b, ^ur) in Eq. ( p0|) is contracted 
to zero, i.e., evolution is "frozen". The renormalization scale should be at least equal 
to the factorization scale, so that the effective coupling has arguments in the range 
controlled by perturbation theory. 

4. Cif(x,y)Q > otherwise evolution to that scale is "frozen" because this re- 
gion is appropriately accounted for by the Sudakov contribution. This helps avoiding 
double counting terms which belong to the resummed rather than to the fixed-order 
perturbation theory. 

5. C i f(x,y)Q > Ci/b; otherwise the two scales /ir = /ip = C%/b are identified in the 
function fuv(x,y). If < A QCD , then f\jy(x,y) is set equal to zero. 

6. Ci/b > Aq CD ; otherwise the function fm(x,y) is set equal to zero. 

To illustrate the difference in technology between approaches employing the conventional 



expression for the full Sudakov exponent |18jl0|, on one hand, and our analysis, on the other 



hand, we show exp(— S) graphically in Fig. ^| for three different values of the momentum 



transfer. In contrast to Li and Sterman []10| , the evolutional contribution is not cut-off at 
unity, whenever C 2 £,Q < G\jb. The dotted curve shows the result for Eq. Q2%D without 
this cutoff. One infers from this figure that their suggestion to ignore the enhancement due 
to the anomalous dimension does not apply in our case because the IR-modified Sudakov 
form factor is not so rapidly decreasing as b increases owing to the IR-finiteness of a™. 
Indeed, as Q becomes smaller, exp(— S) remains constant and fixed to unity for increasing b, 
providing enhancement only in the large-6 region before it reaches the kinematic boundary 
C\jb = Aqcd, where it is set to zero. As a result, for small Q-values, like Qi = 2 GeV, 
the enhancement due to the quark anomalous dimension cannot be associated with higher- 
order corrections to Th, since it operates at larger b- values, and for that reason it should be 
taken into account. Only for asymptotically large Q values, when the IR-modified Sudakov 
form factor and the conventional one become indistinguishable, the evolutional enhancement 
becomes a small effect, strictly confined in the small-6 region, and can be safely ignored. 
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Figure 2. Behavior of the Sudakov form factor with respect to the transverse separation b for three 
different values of the momentum transfer Q\ = 2 GeV, Q2 = 5GeV, and Q3 = 10 GeV with all £j = 1/2, 
and where we set C\ = 2e _7E , C2 = e _1//2 and Aqcd = 0.242 GeV. The dotted curve shows the result 
obtained with af s and Aqcd = 0.2 GeV for Q 2 = 5 GeV using the same set of Cj as before. Notice that 
in this case, evolution is limited by the (renormalization) scale /^r, = t = {m&x^fxy Q, Ci/b}, proposed in 
pCI. However, the enhancement at small b- values due to the quark anomalous dimension is not neglected. 
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Table 1. Different sets of coefficients Cj, and values of the K- factor and the quantity (cf. 
Eq. (p3D) k = In (Cf e 27B_1 /4C|), corresponding to different factorization and renormalization 
prescriptions. 









Scheme parameters Ci 






Choice 




Ci 


C = 


c 3 






canonical 


2exp (-te) 


1/V2 2exp(- 7E ) 


4.565 


-0.307 


SSK § 


exp 


-H 2 7E"1) 


l/y/2 exp 


-|(2 7 E-1) 


2.827 





this work 


2exp (-te) 


exp (-1/2) 2exp(- 7E ) 


exp (-1/2) 4.565 






This behavior of the Sudakov form factor, we emphasize, shows that power-induced 
subleading logarithmic corrections are relevant in the range of currently probed momentum- 
transfer values. Then the advantage of employing such a scheme to calculate observables 
is that the hard contribution to the pion form factor is enhanced. This is because the 
range in which soft gluons build up the Sudakov form factor is enlarged and inhibition of 
bremsstrahlung sets in at larger Q 2 . Let us mention in this context that power corrections 
could also lead to additional suppression of soft gluon emission at large transverse distance 
b. Indeed, Akhoury, Sincovics, and Sotiropoulos j|6| have resummed nonperturbative power 
corrections in b to the hadronic (meson) wave function and found Sudakov-type suppression 
on top of the Sudakov suppression discussed so far. The discussion of such IR-renormalon- 
based contributions in conjunction with our approach would be interesting, but goes beyond 
the scope of the present work. 

Using the techniques discussed above, we obtain the theoretical predictions shown in Fig. 
|3|. Recall that the asymptotic pion distribution amplitude as = 6x(l — x) is employed with 
(3i = 0.883 [GeV~ 2 ] which corresponds to (k 2 ± ) 1/2 = 0.35 GeV (see, 0). A set of constants 
Ci, (i = 1,2,3) which eliminates artifacts of dimensional regularization, while practically 
preserving the matching between the resummed and the fixed-order calculation, is given 
111 table in comparison with other common choices of these constants. In addition, we 
set f(x,y) = ^ s /xy, a choice for the renormalization scale which eliminates large kinematical 
corrections due to soft gluon emission. However, these favored values of the scheme constants 
do not constitute a strict constraint on the validity of the numerical analysis. They merely 
indicate the anticipated appropriate choice of the factorization and renormalization scales 
with respect to observables and theoretical self-consistency. 

Before we proceed with the discussion of these results, let us first present the theoretical 
prediction for the pion-photon transition form factor F 1T1 * 1 {Q 2 ,q 2 = 0) in which one of the 
photons is highly off-shell and the other one is close to its mass-shell. In leading perturbative 
order this form factor is given by the expression (cf. ||26|| ) 



x {AttK (yE bQ) ) e~ s ^ , (35) 
where the Sudakov exponent, including evolution, has the form 

SVy (a:, x, b, Q, C u C 2 , C 4 ) = s(x, b, Q) + s(x, b, Q) - 4r (in) ■ (36) 
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Figure 3. Spacelike pion form factor calculated with (f)^ within our theoretical scheme in comparison 
with existing experimental data |47],^8| . The dotted line shows the enhanced hard contribution to the form 
factor, obtained in our analysis at leading perturbative order and with inclusion of the NLO correction to 
Th (lower solid line). The upper solid line represents the sum of the lower solid line and the dashed line, the 
latter curve giving the result for the soft contribution, computed in ^6|. The dot-dashed line reproduces 
the recent calculation of Tung and Li [T^ ] which does not include an intrinsic fcj_-dependence and uses for 
a s the conventional 1-loop expression. 
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Figure 4. Pion-photon transition form factor calculated with <^> as and IR enhancement (lower solid line). 
The horizontal line represents the asymptotic behavior. The data are taken from B9,M. 



The main difference relative to the previous case is that this form factor contains only one 
pion wave function, whereas the associated hard-scattering part, being purely electromag- 
netic at this order, does not depend on a s . The only dependence on the strong coupling 
constant enters through the anomalous dimensions in the Sudakov form factor. The result 
of this calculation is displayed in Fig. |j. 

All constraints on kinematics set forward in the numerical evaluation of the electromag- 
netic pion form factor are relevant to this case too, except the requirement which deals specif- 
ically with the choice of the renormalization scale, which now is set equal to /ir = C^xQ. 
Another reasonable choice would be //r, = C^y/xEQ, which entails evolution to a lower scale, 
hence reducing evolutional enhancement through r (Ci/b, //r.) by approximately 6%. 

It is obvious from Fig. that the IR-enhanced hard contribution to F V (Q 2 ) is providing a 
sizeable fraction of the magnitude of the form factor. This behavior is IR stable and extends 
from low to high Q 2 values, exhibiting exact scaling. In contrast to the dot-dashed line, which 
is calculated with the conventional ozf\ enhancement is not resulting from approaching 
the unphysical Landau singularity at Q 2 = Aq CD - still appreciable at Q 2 values around 
(6 — 8) GeV 2 - but is induced by the nonperturbative power correction. 

To make these findings more transparent, we compare our results in table ^ with those 
obtained in other analyses. 

Comparison of our values with those calculated by Jakob and Kroll [|llj shows that the 
suppression of the hard part of the form factor due to the inclusion of transverse degrees of 
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Table 2. Calculated pion form factor at two values of Q 2 . The first two columns show the 
results obtained in the present work, in comparison with those calculated by Jakob and Kroll (JK) 
|yj (third column), and by Melic et al. (MNP) @ (last two columns). 



Q 2 [GeV 2 ] 


LO 


LO+NLO 


LO 


LO 


LO+NLO 




(this work) 


(this work) 


(JK) 


(MNP) 


(MNP) 


4 


0.119 


0.167 


0.08 


0.131 


0.211 


10 


0.128 


0.168 


0.08 


0.109 


0.164 



freedom is counteracted by the power- induced enhancement. On the other hand, comparison 
with the values computed by Melic et al. at leading order, by ignoring completely trans- 
verse degrees of freedom, reveals that in the Q 2 domain, where the influence of the Landau 
singularity has died out, there is enhancement. Furthermore, comparing our results with 
theirs at next-to-leading order, we conclude that our choice of the renormalization scale is 
consistent with a proper matching between gluon corrections, calculated on a term-by-term 
perturbation expansion, and those due to the resummed perturbative series. Moreover, the 
scaling behavior of the calculated perturbative contribution to the full pion form factor is 
also improved. 

The calculation of F nl (Q 2 ) provides an additional confirmation of our method. Fig. |] 
shows that our theoretical prediction for this form factor reproduces the recent high-precision 
CLEO im and the earlier CELLO Q] data with the same numerical accuracy as the dipole 
interpolation formula, without using any phenomenological adjustment. 

In this context, we mention, however, that other authors |]5T| , |52| , |53| obtain similarly good 
numerical agreement of F nl * with the experimental data, following different premises on the 
basis of QCD sum rules. 

Another important point worth to be mentioned is that the proposed scheme exceeds 
the quality of the phenomenological predictions derived from schemes which involve satu- 
ration prescriptions [[54]], even if these employ "comensurate rescaling" to "pre-sum" into 



the running coupling constant NLO contributions [p5|] . A more detailed study of the BLM 
procedure []56[| , which helps avoid scheme ambiguities, in combination with our approach 
will be presented elsewhere. 



6. Summary and conclusions 

We have developed in detail a theoretical framework which self-consistently incorporates 
effects resulting from a modification of the running a s by a nonperturbative power correction 
which provides IR universality. Though a deep physical understanding of such contribu- 
tions is still lacking, we have given, as a matter of practice, quantitative evidence that in 
this way it is possible to get a hard contribution to the electromagnetic form factor F n (Q 2 ), 
which is IR-enhanced relative to conventional approaches, using solely the asymptotic form 
of the pion distribution amplitude, hence avoiding end-point concentrated distribution am- 
plitudes. The presented IR-finite factorization and renormalization scheme makes it possible 
to take into account transverse degrees of freedom both in the pion wave function |L1| as well 



as in the form of Sudakov effects [10], without entailing suppression of the form-factor mag- 
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nitude. In addition, use of this modified form of a s (Q 2 ) renders the theoretical predictions 
insensitive to its variation with Q 2 at small momentum values. An appropriate choice of 
the factorization and renormalization scales ensures matching of fixed-order with resummed 
perturbation theory and helps avoid double counting of higher-order corrections due to gluon 
radiation, ensuring this way the self-consistency of the whole perturbative treatment. The 
same procedure applied to F„.o 7 * 7 yields a prediction which agrees with the experimental 
data very well, hence supporting the theoretical basis of our approach. 

Taking our results for the pion form factor at face value, they seem to suggest that its 
magnitude may be significantly larger than predicted in previous analyses, even without 
employing broad pion distribution amplitudes. However, because of the restricted accuracy 
of the existing experimental data such a conclusion would be surely premature. 

We believe that the insight gained through our analysis gives a strong argument that 
if power corrections to the effective coupling of QCD really exist, their implications and 
applications are important and revealing. In this respect, the Shirkov-Solovtsov approach 
may provide a convenient tool for improving theoretical predictions for exclusive observables, 
based on perturbation theory. 
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